NUMERICAL MODELLING OF CORRELATIONS 
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We demonstrate algorithm for numerical modelling of Bose-Einstein correlations 
(BEC) formulated on quantum statistical level for a single event and exploring the 
property that identical particles subjected to Bose statistics do bunch themselves 
in a maximal possible way in the same cells in phase-space. 



Bose-Einstein correlations (BEC) between identical bosons are supposed 
to provideinformation on space-time development of multiparticle production 
' ijrocessesEl. The usual Monte Carlo event generators modelling such processes 

'nI ^ 0^ because their probabilistic structure, excludes a priori the genuine BEC, 

which are of purely quantum statistical origin^ JZLne can only model BEC by (i) 
suitably changing output of these generators la crEl or by (m) building generator, 
which properly incorporates the bosonic character of produced particles u. In 
both cases the goal is to reproduce the experimental two-particle correlation 
p • function C2[Q = \pi - Pj\) = N2{pi,Pj)/Ni{pi)Ni{pj). In (i) this is achieved 

by suitable bunching of the finally produced identical particles in phase-space 
performed using special weights constructed from the output of the event 
, generator. In (m) the particles are already being produced in properly bunched 

way by means of special generator constructed using specific statistical model 
(providing Bose-Einstein or geometrical distribution of particles in each bunch 
■ which is identified with a single emitting cell in phase space Whereas (i) 

^ I can be applied only for all events and is (via weights) sensitive to the space- 

time structure of the production process provided by event generator, the 
(m) applies already on a single event level but its generator bears no a priori 
information on the space-time structure of the production process, it uses 
instead nonstatistical character of fluctuations it produces. 

We propose generalization of the second approach to make it applicable 
also to other generator^. To better understand our reasoning let us remind 



"^Cf. Q for more details, especially in what concerns the hadronization model CAS used in 
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that classical weight method □ amounts to multiplying each event by spe- 
cial weight, i.e., event is counted many times if it already possesses, by shear 
chance, traces of desired bunching. In terms of philosophy of which we shall 
follow, it could be seen as selecting events already possesing (to some degree, 
at least) a more bosonic character than other ^events (i.e., in which particles 
are already bunched in way resembling that of a) . What we propose is similar 
approach but performed already on a single event level. Namely we propose 
to search for the bosonic configurations of particles existing already in each 
event because of the internal nonstatistical fluctuations provided by event 
generator. Namely, there are groups of particles resembling those obtained in 
U, modulo only the fact that they usually have different charges allocated to 
them whereas particles in □ are of the same charge. We propose therefore to 
endow such bunches of particles with the same charge to an extent limited only 
by the overall charge conservation. This means that in cases where charge al- 
location has been already provided by event generator we shall neglect it and 
perform new charge allocation keeping, however, the total number of particles 
of each charge the same as given by this generator. Notice that we do it for 
each single event, keeping intact both the original energy-momentum pattern 
provided by event generator (i.e., conserve the energy-momentum) and all in- 
clusive single particle distributions. Leaving those interested in more details 
to l3 we shall only say that to get desired result it is enough to select one of 
the produced particles, allocate to it some charge, and then allocate (in some 
prescribe way) the same charge to as many particles located near it in the 
phase space as possible. In this way one forms a cell in phase-space, which 
is occupied by particles of the same charge only. This process should then 
be repeated until all particles are used and it should be such that one gets 
geometrical (Bose-Einstein) distribution of particles in a given cell. This pro- 
cedure changes the charge flow pattern provided by event generator^ retaining, 
however, both the initial charge of the system and its total multiplicity distri- 
bution. The procedure of formation of such cells is controlled by probability 
P that given neighbor of the initially selected particle should be counted as 
another member of the newly created emmiting cell in phase spac^. 

Referring to Q for more details we shall illustrate in Figs, la and lb 



calculating results in Fig. 1. 

''It amounts to allowing formation of multi-likc-charged object on intermediate steps of 
hadronization process. Therefore this method works only when such possibility exists in a 
given generator. 

'^It is important to realize that, because we do not restrict a priori the number of particles 
which can be put in a given cell, we are automatically getting BEC of all orders (even if we 
use only two particle checking procedure at a given step in our algorithm). It means that 
C2{Q = 0) calculated in such environment of the possible multiparticle BEC can exceed 2. 
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our attempts to describe (separately) e+e" dataJpy DELPHI on BEC □ and 
intermittency U using the so called CAS model Il3 (see also Q) whereas Figs. 
Ic and Id show the respective intermittency and BEC obtained when using 
parameters from the fits above. The results, although not totally satisfactory, 
are encouraging given the simplicity of CAS model used. Application of our 
method with other, more sophisticated event generators should answer the 
question of its final applicability. 
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Figure 1. The examples of best fits to data on e+e~ annihilation by DELPHI on BEC * 
in (o) and, separately, on factorial moments ^ Fi in (6) (M is number of bins) using simple 
cascade hadronization model CAS (cf. '^'^^ for details). In (c) are shown factorial moments 
for parameters used in (a) (when fitting BEC) whereas in (d) BEC for parameters used in 
(6) (when fitting factorial moments). Notice that left panels ((o) and (c)) are obtained for 
3 sources whereas right panels ((6) and [d)) for 2 sources. To fit Fs one needs P = 0.75 
and two sources, but in this case the calculated F2 overshoots data by ca 50%. 
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